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Most studies of PT-symmetric quantum- mechanical Hamiltonians have considered 
the Schrodinger eigenvalue problem on an infinite domain. This paper examines the 
consequences of imposing the boundary conditions on a finite domain. As is the 
case with regular Hermitian Sturm-Liouville problems, the eigenvalues of the VT- 
symmetric Sturm-Liouville problem grow like n 2 for large n. However, the novelty 
is that a VT eigenvalue problem on a finite domain typically exhibits a sequence 
of critical points at which pairs of eigenvalues cease to be real and become complex 
conjugates of one another. For the potentials considered here this sequence of critical 
points is associated with a turning point on the imaginary axis in the complex plane. 
WKB analysis is used to calculate the asymptotic behaviors of the real eigenvalues 
and the locations of the critical points. The method turns out to be surprisingly 
accurate even at low energies. 
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I. INTRODUCTION 

This paper uses WKB analysis to examine the approximate solutions of complex non- 
Hermitian "PT-symmetric Sturm-Liouville eigenvalue problems on finite domains. These 
problems are qualitatively different from ordinary Hermitian eigenvalue problems because, 
as we will show, there is a sequence of critical points at which the eigenvalues become 
pairwise complex. WKB provides an extremely accurate asymptotic approximation to the 
locations of these critical points. 

A conventional Sturm-Liouville eigenvalue problem in Schrodinger form is a second-order 
differential equation 

-ip"(x) + V(x)if)(x) = \ip(x), (1) 
where A is the eigenvalue, accompanied by a set of homogeneous boundary conditions 

V>(a) = 0, i/>(b) = 0. (2) 

If a and b are finite and the potential V(x) is real and smooth for a < x < b, then this 
eigenvalue problem is said to be a regular Sturm-Liouville problem. 

WKB theory gives a good approximation to the solution of this problem for large eigen- 
values A pQ. Equation ^ takes the form ip"(x) + Q(x)ip(x) = 0, where Q(x) = A — V(x), 
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and if A ^ 1, we may assume that Q(x) 7^ on the interval a < x < b. Thus, there are no 
turning points. When A> 1, the WKB approximation satisfying ip{a) = 0, 



^ xH [gRF sin 



ds ^\ - V(s) 



(3) 



is uniformly valid over the entire interval. Imposing the boundary condition ip(b) = then 
gives the eigenvalue condition 

/ ds y/X n - V(s) ~ mr (n= 1,2,3,...). (4) 

J a 

For large A this gives an accurate asymptotic approximation to the eigenvalues: 

T7?7T 2 

A„~^-^ (n->oo). (5) 

Note that the eigenvalues grow like n 2 for large n. This result is general and holds for all 
regular Sturm-Liouville eigenvalue problems because to leading order in the WKB approx- 
imation the eigenvalues become insensitive to the potential V(x) and thus the eigenvalues 
approach those of a square-well potential. 

In this paper we study the complex "PT-symmetric version of the eigenvalue problem in 
([!]) and ([2]). Now, instead of the potential being real, the eigenvalue problem takes the form 

— ip"(x) — gV(ix)ip(x) = \ip(x), (6) 

where g is a coupling constant and V(x) is a real function of its argument. In order to 
respect the VT symmetry, the boundary conditions are imposed on the real-x axis at parity- 
symmetric points: 

^(±1) = 0. (7) 

We will see that WKB provides an excellent asymptotic solution to this problem for large A 
but that one-turning-point analysis is required. 

Eigenvalue problems like this have already been studied numerically in the literature in a 
variety of physical contexts. An example of a PT-symmetric Hamiltonian having a complex 
Sturm-Liouville eigenvalue problem like that in (J6]) and ^ was first discussed by Giinther, 
Znojil, and Wu [2] in the context of the Squire equation of hydrodynamics. Essentially 
the same Hamiltonian and boundary conditions occur in the context of superconducting 
wires [3], and again in a different guise in the consideration of the magnetic resonance signal 
of spin-polarized Rb atoms near the surfaces of coated cells 0j. The Schrodinger eigenvalue 
equation of Ref . [3] , 

- if/'(x) - ilxtp{x) = AV>0), (8) 

is posed on the finite domain |x| < 1 with homogeneous boundary conditions ^(±1) = 0. 
The eigenvalues A are plotted in Fig. [T]as functions of the real coupling constant /. Observe 
that when < / < 12.31, the eigenvalues are all real. This parametric region is called the 
region of unbroken VT symmetry. However, as / reaches the critical value 12.31, the two 
lowest eigenvalues become degenerate, and as / increases past 12.31, these eigenvalues split 
into a complex-conjugate pair. Thus, for / > 12.31 the eigenspectrum is no longer entirely 
real. As I continues to increase, more and more pairs of real eigenvalues become degenerate 
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FIG. 1: (Color on line) Real parts of eigenvalues of the 'PT-symmetric eigenvalue problem 
taken from Ref. [3]. 



and split into complex-conjugate pairs. The critical values of a coupling constant at which 
the eigenvalues become degenerate are often called exceptional points [5]-|10j. 

Another example of a PT-symmetric Hamiltonian giving rise to a Sturm-Liouville eigen- 
value problem on a finite domain is 



H 



de 



2 +ig cos 9, 



(9) 



where g is a real parameter [TTj. Here, the eigenfunctions are required to be 2n periodic 
and odd in 9. The region of unbroken VT symmetry is \g\ < 3.4645 pj] (see Fig. [2]). Note 
that the behavior of the eigenvalues is virtually identical to that shown in Fig. [IJ It is this 
feature that motivated us to use WKB analysis to investigate these problems in a more 
general context. 

This paper is organized simply. In Sec. [IT] we describe the WKB calculation, culminating 
in Eq. (32), which is our principal result. Sec. Ill gives numerical results based on (32) for 
a variety of Hamiltonians, including in particular those of ^ and ([9]). Section IV contains 
some brief concluding remarks. 
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FIG. 2: Real and imaginary parts of eigenvalues of the PT-symmetric Hamiltonian ([9]), taken from 
Ref. The region of unbroken VT symmetry is \g\ < 3.4645. 

II. WKB CALCULATION OF EIGENVALUES AND CRITICAL POINTS 

Our objective is to find accurate asymptotic approximations to the eigenvalues A in ^ 
for large A and hence to determine approximately the critical values of g. In our WKB 
approximation we treat both the eigenvalues A and the associated values of g as large and 
proportional. Thus, we begin by substituting 

A = ag. (10) 

The eigenvalue equation ^ then takes the form 

r(x) = -gQ(x)iP(x), (11) 

where Q(x) = V(ix) + a. 

We will see that for large A, the asymptotic solution to (11) is controlled by a turning 



point on the imaginary- x axis at x = ib, where b satisfies the equation 

V(-b) + a = 0. (12) 

Near the turning point at ib, we write x as x = ib + yc, where ?/ < 1. Thus, a one-term 
Taylor approximation to Q(x) near ib is Q(x) ~ iycV'(-b), and the differential equation 



(11) is approximately 

jz^y) = ^gc 3 V'(-b)y4>(y). (13) 

We then convert (13) to the standard Airy differential equation 

d 2 

-j^i>(y) = y*P(y) (14) 

by setting 

c = 7 e m/6 , (15) 

where 7 = [gV'(— 6)] -1 / 3 is real. By choosing the sign of g appropriately, we may take 7 to 
be positive. The relationship between the x and y variables is given explicitly by 

y = ^L^ e --/6. (16) 
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Two linearly independent solutions to the Airy equation (14) are Ai(y) and Ai(uy), where 
to = e 2m / 3 is a cube root of unity. Hence, the general solution to ( 14 ) is 

i/>(y)=K 1 M(y) + K 2 M(uy), (17) 

where K\ and K% are arbitrary constants. Thus, near the turning point at x = ib on the 
imaginary axis, the solution ip(x) to the Schrodinger equation (11) is 

x — ib 



- A",Ai ( ^V"^ 6 

7 



+ A" 2 Ai u 



in/ 6 
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(18) 



Away from the turning point at x = ib the solution to (11 ) can be expressed in terms of 
WKB functions because g ^> 1. When Rex < 0, the WKB approximation to ip{x) on the 
left takes the form 



\Q{ x )] l/ ' 



exp 



ib 



ds y/gQ(s) 



+ 



[Q{x)]W 



exp 



— i 



and when Rex > 0, the WKB approximation on the right is 



Ri 



[Q(x)]V4 



exp 



ds \/gQ(s) 



ib 



+ 



R-2 



exp 



ds y/g~Q(s) 



ds \/gQ(s) 



it, 



(19) 



(20) 



where L\, L 2 , R\, and R 2 are arbitrary constants. Note that the sense of integration in (19) 
and (20) is from left to right; in (19) the integration is rightward and towards the turning 
point at ib and in (20) the integration is rightward and away from the turning point. 

We must now find equations that relate the six constants in the approximations to ip{x) 
in (17) - (20). Imposing the boundary condition ^(—1) = on ^l(x) in (19), we obtain a 

ds \/gQ(s) 



condition relating L\ and L 2 , 



= Li exp 



ib 



ds \/g~Q{s) 



-i 



+ L 2 exp 



— i 



ib 



-i 



(21) 



and imposing the boundary condition = on iPr(x) in (20), we obtain a condition 

relating Ri and R 2 , 



= Ri exp 



ds y/ gQ(s) 



ib 



+ R 2 exp 



—i 



ds y/gQ(s) 



ib 



(22) 



Four additional conditions relating the constants can be found by matching asymptoti- 
cally the WKB approximations (19) and (20) to the Airy approximation (17) or (18) near 
the turning point. To perform the match, we must show that in an overlap region near the 
turning point, further asymptotic approximations to each of the asymptotic approximations 
that we have already found are identical. We will need to use two approximations to the 
Airy function for large argument that are valid in the appropriate Stokes' wedges PQ: 



Ai(z) 



1 



20F 



and 



Ai(-z) 



~ —j=Z 



-1/4 



exp I --z 



2 3/2 , ^ 



3/2 



(\z\ > 1, |arg z\ < 7r) 



sin I -z 



\z\ 3> 1, lart 



-z < 



2tt 



(23) 



(24) 
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We will perform the asymptotic match in the y variable. Because y = (x — ib)/c and 
c = O (g -1 / 3 ) is small for large g, it follows that when x is near ib, y may be treated as large. 
Thus, it is valid to use the asymptotic approximations (23) and (24) for the Airy functions. 



First, we examine the WKB approximation i(jr(x) in (20). This WKB approximation is 
valid in the right-half x plane away from the turning point at ib; that is, outside the small 
circle in Fig. [3} From the asymptotic approximations 



ds \/ gQ{s) 



-1/4 



ff l/41/2 e -i»/6 1/ -l/4 ) 



(25) 



the asymptotic approximation to the WKB approximation to ipn(x) in (20) becomes 



1/4 1/2 -M-/6 -1/4 



V 



Ri exp 



3/2 



+ # 2 exp ' -»< 3/2 



(26) 



This approximation is valid as we approach the circle from the outside. 




FIG. 3: Schematic path for the WKB integral from the end points x = ±1 via the turning point 
at x = ib. Inside the circle, very close to the turning point, the WKB approximation must be 
matched on either side with the asymptotic approximation to the Airy solution. 



The left panel of Fig. [4] shows the complex-x plane and the right panel shows the cor- 
responding complex-?/ plane in the vicinity of the turning point at x = ib. From (16), we 
see that the complex-?/ plane is rotated by —30° relative to the complex-x plane. The Airy 



functions in (17) are oscillatory along the wiggly lines in the right panel. 



The Airy function approximation in (17) is valid inside the circle in Fig. pi We must 



match the WKB approximation in ( 26 ) to the Airy approximation along the solid line (R) 
in the left panel of Fig. |1J This line corresponds to the two wiggly lines marked (R) in the 
right panel. Because both wiggly lines satisfy the conditions of (23), as we approach the 



right edge of the circle from inside, we use only this asymptotic approximation to obtain 



ip(x) 



-1/4 



71 " 



A>xp( --y 3/2 



(27) 
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FIG. 4: Matching paths in the complex- a; plane (left panel) and y plane (right panel). 



This asymptotic match produces two algebraic equations for the coefficients: 



20T 

20F 



(28) 



Next, we further approximate the WKB approximation in (|19|). Again, using the formulas 

2 



in (25), we find that in the left-half x plane 

Li exp 



^l/4 T l/2 e i W /3(_ 2/ )-l/4 



-yf' 2 



f L 2 exp ( --i{-yf 2 



(29) 



For the Airy approximation (17) we must now perform the asymptotic match along the solid 



line (L) in the left panel of Fig. |4j This line corresponds to the two wiggly lines marked (L) 
in the right panel. In this case the wiggly line along the negative-?/ axis requires that we use 
the asymptotic approximation (24) for the Airy function Ai(y). For the other Airy function 
Ai(uy) we use the asymptotic approximation (23). This allows us to obtain the following 

(30) 



asymptotic approximation to the Airy-function approximation to i/j(x) in (17): 



71 



i-y) 



-1/4 



sin 



{-yf 12 + 



71 



4 



+ _^(_ y )-l/4 e Vl3 



2Jn 



exp 



-y) 



3/2 



This asymptotic match produces two further algebraic equations 



^2^ 1/4 7 1/2 e- w/6 , 



2Jn 



(31) 



Finally, we combine all six algebraic equations (21), (22), (28), and (31), and obtain a 
secular equation that determines the eigenvalues: 



sin 



ds y/gQ(s) 



+ - exp 



ib 



ds \/gQ(s) 



ds y/gQ(s) 



ib 



0. 



(32) 



8 



Equation (32) is our main result. This equation is real, which is a general feature of all 
PT-symmetric secular equations [12] . To see that it is indeed real note that the argument 
of the sine can be written as 2ReJ + , while that of the exponential can be written as 2 Im/ + , 
where J+ = f. b ds ^ 'gQ(s) 

The first term of (32) is what one would obtain from a path going directly along the 
real axis from x = — 1 to x = 1 without going through the turning point at x = ib. This 
corresponds to the no-turning-point result of ^ for the Hermitian case. We shall see in 
the next section that this term by itself gives very little structure. The second term is the 
exponential of a generically large real number. When that number is large and negative it 
makes very little difference to the calculation of the eigenvalue, and when it is large and 
positive the equation has no real solutions because | sin#| < 1. All the interesting structure, 
including the critical points, comes from the interplay of the two terms in the region where 
the exponent passes through a zero. 



III. NUMERICAL CALCULATIONS 



In Fig. |5| we show the predictions of (|32|) for the real eigenvalues of the Airy potential 

V{ix) 



ix occurring in (|8]) (with g taking the role of I). These are the solid lines, which 
are essentially indistinguishable from the numerical results of Fig. [TJ Remarkably, the range 
of validity of the WKB approximation extends to small values of g and A. In Fig. [5] we also 
show as dotted lines the result of the no-turning-point WKB approximation, namely the first 



term of (32). While this approximation reproduces correctly the square- well eigenvalues for 
0, it fails to reproduce the interesting structure, which arises as a result of the interplay 
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between both terms of (32). For this potential the full WKB approximation predicts that 

g/y/3, which is actually an exact result [2| IT3~ 



the critical points occur at A 

In Fig. [6] we show the analogous results for the sinusoidal potential in (9|) (with 9 
7i / 2 — x). The same features of the approximation are true here too, and (32) reproduces 



veryclosely the numerical results of Fig. [2j The reason that the spectra of the two potentials 
in (81) and (9) are so close is that within the range of the integrals J ds >JgQ(s) occurring 
in (32) the function sins is well approximated by s. Note that in this case b 



arcsm a 



is multivalued, so there is a series of turning points on the imaginary axis. For the WKB 
calculation we considered only the nearest turning point, which is close to ia. 

In Fig. [7] we show the analogous results for the potential V(ix) = isin(2x). Again, these 
are indistinguishable from numerical results obtained from a shooting method, but now they 
differ from the eigenvalues of the scaled Airy potential V(ix) = 2ix, which demonstrates 
that within the range of the integrals J ds^gQ(s) it is not valid to neglect the s 3 term in 
the expansion of sin(2s). 

What happens if we apply our WKB approximation to the potential V(ix) = — ix 3 ? Now 
there are three complex turning points, which are located at x — ib{l,u,u> 2 }. Since our 
WKB analysis involved only one turning point, we chose a path going through the turning 
point at ib. The numerical results in Fig. M show a single low-energy critical point, which 



the WKB approximation fails to reproduce. In this case the second term of (32) is always 



negligible, so there is no effective interplay between the two terms. On the other hand, the 
first term tracks very well the curves for the higher energy levels. It is conceivable that 
a WKB path going through the other two turning points would reproduce the low-energy 
structure, and we intend to address this problem in a future publication. 
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FIG. 5: WKB approximation for the real energy levels of the Airy potential V(ix) = ix. The 
no-turning-point approximation (dashed line), which is obtained by neglecting the second term 



of (32), exhibits no critical points. The full result (solid line), which includes the second term, 



indeed gives the critical points as a consequence of the interplay of the two terms. The solid line 
is essentially indistinguishable from the numerical results in Fig. [TJ even at low energies. 
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FIG. 6: WKB approximation for the real energy levels of the potential V(ix) = is'mx. The 
notation is as in Fig. [5] The solid line is essentially indistinguishable from the numerical results in 
Fig. [2j even at low energies. 




FIG. 7: WKB approximation for the real energy levels of the potential V(ix) = isin(2x). These 
are distinct from those of the scaled Airy potential V(ix) = 2ix. 



X 

100r 




FIG. 8: (Color on line) WKB approximation for the real energy levels of the potential V(ix) = —ix 3 
compared with the numerical results (red dots). The numerical results show just one low-energy 
critical point in the range of g considered, while the WKB approximation gives no critical points. 
This is because the second term of (32) is negligible in this case. 
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IV. COMMENTS AND DISCUSSION 



In this paper we have used WKB analysis to derive the simple formula (32), which 
gives an extremely accurate approximation for the energy levels A of PT-symmetric Sturm- 
Liouville eigenvalue problems. This derivation requires the use of one-turning-point analysis. 
If a no-turning-point analysis is used one still obtains a good approximation to the high- 
lying energy levels for fixed coupling constant g. However, the no-turning-point analysis is 
unable to reproduce the critical points that occur when g and A are both large. What is 
most remarkable is that the one-turning-point formula gives an accurate description of the 
spectrum and critical points even when g and A are not large. 

Particular examples that we have studied included the ix and i sin x potentials. In these 
cases we were able to understand the close equality of their respective spectra and why the 
spectra of the pair of potentials 2ix and i sin(2x) were not equal. In the case of the ix 3 
potential the WKB approximation correctly reproduces that higher energy levels, which do 
not exhibit any critical behavior. There is, however, a low-energy critical point, as seen in 
Fig. |8j which our WKB approximation does not reproduce. It is an open question whether 
a two-turning-point approximation would be able to do so. 
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